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Abstract 

The antisymmetrised two-electron functions are constructed for the point group T^, i.e. the 
symmetry group of unconventional superconductor PrOs4St>i4., and its subgroup D2^.The nodal 
structure of these function depending on the wavevector group is investigated. Theoretical nodal 
structure of these functions made possible to explain the experimental nodal structure PrOs4Sbi4. 
as a result of time-reversal symmetry violation. 
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I. INTRODUCTION 



The superconducting transition is accompanied by different types of symmetry violations. 
In particular these symmetry violations appears as anisotropic structure of the supercon- 
ducting order parameter ( SOP), vanishes at lines and points of Fermi surface [1]. The SOP 
in many of unconventional superconductors as antisymmetrical (UTh)Bei3 [2], S^RuO^ 
[3,4] and PrOs^Sb^ [5-9] is triplet whose spatial part is odd. The muon spin relaxation 
experiments in the three above material magnetic fields, which result in time-reversal sym- 
metry violation[10]. 

Thermal conductivity and specific heat measurements in magnetic field, indicated that 
SOP of FrOs 4 S'6i2depends on the magnitude of the field [11] and it was concluded that 
at high fields SOP has six point nodes in the directions [±1,0,0], [0, ±1,0] and [0,0, ±1] 
and at low fields SOP has two point nodes in the directions [0,0, ±1]. Hence it follows 
that PrOs^Sbi2 is the first heavy-fermion superconductor having two phases with point 
nodes and a transition between them occurs in superconducting state. It was suggested 
that superconductivity in PrOs^Sb^ is connected with electric quadrupole interaction in 
distinction to other heavy-fermion superconductors. 

The nodal structure of SOP is closely connected with interactions responsible for Cooper 
pairing e.g. magnetic interactions instead of electron-phonon interactions. Experimental 
investigations of SOP nodal structure and theoretical interpretations of these experiments 
are basic to an understanding of the nature of such interactions. There are three group 
theoretical approaches to the SOP investigations based on point groups, continuous groups 
and space groups. The point-group approach to SOP [1,12,13] is based on the reduction 
of spherical functions, corresponding to isotropic case onto the crystal point group. The 
choice of basis functions of irreducible representation (IR) of point group is not unique 
and hence the conclusions of this approach on the nodal structure of SOP are ambiguous 
[14] . According to the results of point group approach there are no symmetry reasons for 
vanishing of triplet SOP on the lines of Fermi surface [13] (Blount theorem). This statement 
is connected with additional degree of freedom of spin part of the wavefunction. In this 
situation the majority of triplet superconductors, in which SOP vanishes on lines of Fermi 
surface , belong to the exclusion of this theorem. We will show below that in order to 
explain these exclusions, the symmetry reduction of spin part , connected with time-reversal 
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symmetry violations should be taken into account. 

The continuous group approach unites singlet SOP and antiferromagnetic order parame- 
ter by the operators of SO (5) group [15]. In this approach a two-fermion state is described 
by the five dimensional superspin with following components: singlet order parameter for 
particles, singlet order parameter for holes and three components of antiferromagnetic order 
parameter. 

The space-group approach to SOP [16-23] is based on the Ginzburg-Landau [24] statement 
that the symmetries of SOP and of the Cooper pair wavefunction are identical, on the 
Anderson symmetry description of singlet and triplet pairs in a general point in k- space 
and on the induced representation method [26,27]. In the present work making use of the 
space-group approach general triplet two electron wave-functions are constructed for group 
Th and its subgroup D 2 h- The latter group corresponds to the symmetry of triplet spin in 
the case of time-reversal symmetry vilation. This resultes are used to identify the symmerty 
of phases with different nodal structure of PrOs^Sb^. 

II. METHOD OF CALCULATIONS AND RESULTS 

In spherically symmetric case two electrons with opposite momenta are connected into 
a Cooper pair. According to Anderson [25], the wavefunctions of electrons in singlet pair 
are connected by time reversal 6*,and three components of triplet pair are written in second 
quantization form as: 



where / is the space inversion. In L-S coupling singlet and the three components of triplet 
two-electron functions are written as: 



{c\c\ k ) 



(1) 



(2) 




(3) 



(4) 
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tf* = $X.™= -1,0,1 (5) 

where S° and are singlet and triplet functions. Corresponding spatial parts are given 
by two following formulae: 



(6) 



(7) 



where subscripts 1 and / correspond to an arbitrary chosen wave vector in a Brillouin 
zone and the result of action of space inversion on it, and superscripts denote the number 
of electron coordinate. Owing to translation symmetry one-electron functions in crystal are 



Correct translation invariant two-electron wavefunction in crystal is expressed as antisym- 
metrical linear combination of the wavefunction belonging to all prongs of the star. It follows 
from the induced representation theory, [21,26,27] that these functions may be constructed 
by the action of projection operators on the functions (4) and (5). It follows from the reci- 
procity theorem [27], that the projections are non zero for those IRs of the whole group , 
whose characters are not orthogonal to characters of IRs of the subgroup H + IH, corre- 

— * 

sponding to function (4) and (5). Hence it follows that for k a general point in a Brillouin 
zone all even IRs are possible for singlet pair and all add IRs are possible for triplet pair and 
the number of appearance of each IR equals to its dimension. If the wave vector group// 
consists of more then one element, the characters of IRs for singlet and triplet pairs are 
given by two following formulae: 



where h is an element of H, D is IR of H, the plus sign corresponds to singlet pair and 
minus sign corresponds to triplet pair and S is an element of double coset in decomposition G 
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characterized by the group H of the wave vector k and the star of the wave vector {k}. 




(8) 



X ± (Sh) = X (D(5h5h)) 



(9) 



with respect to H. Since the total momentum of a Cooper pair equals zero, we consider the 
double cosets defined by 5 reversing the direction of k. For k a general point of a Brillouin 
zone 5 is a space inversion only and in symmetrical direction we can also take as S a 180° 
rotation around the axis perpendicular to k. 

For real one-dimensional IRs of H we obtain from formulae (8) and (9) the unit represen- 

— * 

tation of H + SH for singlet pairs. For triplet pairs characters of all elements unchanging k 

— * 

equal to 1 and characters of all elements changing direction of k equal to -1. It is clear that 
in general case this two IRs don't exhaust all IRs of H + 5H. It can be easily show that if 
IRs of H are two-dimensional, there are singlet pairs with odd spatial part and triplet pairs 
with even spatial part [16]. Thus formulae (8) and (9) generalize Anderson approach for 
the case of symmetrical directions and planes in a Brillouin zone. 

Let us consider the case of k on the plane of symmetry, i.e. H = C s . We easily obtain 
that spatial part of singlet pair belongs to IR A g of C 2 /i and that of triplet pair belongs to IR 
B u of C 2 h- In the case of singlet pairs the absence of another even IR B g implies, that some 
of the even IRs of the whole group will be forbidden for two-electron states on this plane. 
For triplet pairs without spin the absence of another even IR also implies that some odd IRs 
will be forbidden at this plane. But the account for the spin lifts this restriction. In the case 
of time-reversal symmetry one can take real triplet spin basis functions, which are usually 
denoted as x, y and z. On the symmetry planes these functions belong to two different 
even IRs of group C 2 h and their products on the IR of spatial part B u include all IRs of 
C 2 h- Hence it follows that there are no symmetry requirement for vanishing of triplet SOP 
on the planes of symmetry (Blount theorem[13]). The existence of internal magnetic field 
results in time-reversal symmetry violation. In this case different directions of spin become 
non-equivalent and it is natural to suppose that only one of them corresponds to Cooper 
pairing. In this case the number of possible IRs on the planes of symmetry is reduced and 
symmetry restrictions fro IRs of the whole group become possible. Let one of IRs, forbidden 
on the plane of symmetry corresponds to Cooper pairing. The intersection of this plane 
with the Fermi surface results line on nodes. If the IR, corresponding to Cooper pairing 
forbidden on the line of symmetry, the intersection of this line results a point node of SOP. 

Comparing theoretical nodal structure of SOP with experimental one, we can identify 
IRs corresponding to Cooper pairing [18,19]. Group theoretical analysis of possible IRs in 
k z direction for T h group is carried out in Table 1. It follows from this table, that for singlet 
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pairs all even one-dimensional IRs, but three dimensional IR is forbidden. For triplet pairs 
without account for spin three dimensional only odd three dimensional IR is possible, but 
all one-dimensional IRs are forbidden. When the spin is included in the whole symmetry, 
all IRs are possible, but the type of IR depends on the orientation of the spin relative to 
Z-axis. Triplet basis functions for the case of k z directions are presented in Table 2. It 
should by noted that despite of similarity of these functions to that of ref. [28], physical 
significance is different. In refs. [1,12,13,28] basis vectors k x , k y and k z change continuously 
and one formula (e.g. formula (12) of ref. [28]) describes all directions of k. On the other 

— * 

hand the data of Table 2 correspond to to the directions of coordinate axes only. For k a 
general point in a Brillouin zone the shape of our basis functions is different and we'll obtain 
it by projection operator technique. 

In general point of a Brillouin zone the dimension of the wave vector star equals to the 
number of point group elements, and the number of basis functions in each case (singlet or 
triplet) is two times less, since the inversion was already used in the construction of functions 
(6) and (7). Triplet basis function for IRs A u and T u constructed by projection operator 
technique are presented in table 3. The nodal structure of basis functions of E\ u and E 2u is 
the same as for A u and these functions are not shown in this table. Function for k\ a general 
point in a Brillouin zone is denoted as §\ and the results of action on it by operations h 2 , 
h 3 and h 4 [29] ( 180° rotations about axis X, Y and Z respectively) are denoted as <3>| 
and $4. The rotations h 5 and /19 transform Z axis into positions Y and X respectively and 
function of last row of IR T u into the second and the first row. There is single spatial basis 
in the cases of IRs A u , E lu and E 2u . When the triplet spin in included the dimension of 
each basis is increasing in three times. In the case of IR T u there are three independent 
spatial basis sets, which can be obtained starting projections from different rows. These 
basis sets are denoted in Table 3 by superscripts a, (3, and 7. It should be also taken into 

— * 

account that in each point in k- space there are three directions of spin are possible. Thus 
in the case of IR T u we have 9 independent basis sets. This dimension follows from the fact 
that total dimension of two electron basis equals to the square of dimension of one-electron 
basis. We show only four of this basis sets T"(z), T£(y), T®(z) and T£(z). If the wave 
vector k\ is approaching one of the coordinate axis four vectors connected by 180° about 
this axis and reflections in planes normal to other axis coordinate axis merge and there are 
two possibilities: basis functions transform into the function presented in Table 2 or cancel. 
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The latter case corresponds to a point node of the wave functions. The nodal directions also 
shown in Table 3. Different basis sets for each of IRs A u and T u vanish at different directions. 
Since in the absence of external magnetic field all basis sets are physically equivalent, there 
are no symmetry requirements for vanishing of pair functions in the symmetry directions. 
This conclusion generalizes the Blount theorem for the case of symmetry directions in which 
the dimension of IR of wave vector group equals to unity. 

The above conclusion implies that experimentally observed point nodes of SOP of 
PrOs^Sbu are connected with symmetry violations. Experimentally observe magnetic fields 
in PrOs^Sbu result in time-reversal symmetry violations [10]. In the presence of magnetic 
fields total spin projections M s —1 and M s =-l are not equivalent to each other and its im- 
possible to construct their real combinations x and y. In this case the triplet two-electron 
state splits into three states, corresponding to tree components of total spin. According 
to Anderson [25] the wavefunction of Cooper pair has a time -reversal and space inversion 
symmetry. Singlet pair and M s =0 pair (in our notations z) component of triplet invariant 
with respect to these operations. Time-reversal interchanges two energetically non equivalent 
states M s —1 and M s —-1 and these state should be excluded from consideration. Also spatial 
operations changing the direction of z component of spin should be excluded. Thus in the 
case of time-reversal symmetry violation z component of triplet pair has a D 2 h symmetry. 

Basis functions of z component of triplet pair in D 2h symmetry are presented in Table 4. 
In a general point of a Brillouin zone all IRs are possible for such a pair, but different IRs 
reveal different structure of point nodes at coordinate directions. Two-electron function of 
Au symmetry has four nodes in directions ±X and ±Y. Functions B lu and B 2u have nodes 
in directions ±X, ±Z and ±F, ±Z respectively. Since in the case under consideration Z is 
the direction of magnetic field, these two functions are physically equivalent and only nodes 
in directions ±Z follow from the symmetry. Finally, function of B 3u symmetry has 6 nodes 
in directions of all coordinate axis. Hence, there are tree possible nodal structures of SOP: 
Au with four nodes, B\ u +B 2u with two nodes and B 3u with six nodes. 

III. DISCUSSION OF THE RESULTS 

The phenomenological model p + h function was proposed in ref. [30] for A-pase in 
PrOs A Sb l2 : 
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A A (k) = A(l -k\-k\- kt) 



(10) 



This formula describes six nodes proposed in ref. [8]. This nodal structure corresponds 
to the B Su pair function of Table 4. For the B— phase with nodes in Z— axis following 
formula was proposed in ref. [30]: 



This nodal structure corresponds to the sum of IRs B\ u and B 2u of Table 4. The model 
of of functions s- and d— symmetry was proposed in ref. [31]. This model describes the bi- 
phase by the anisotropic s- function with six nodes and B- phase, by the function s+id x 2_ y 2 . 
In ref. [32] both states were described by the function s+g type. Also the state with six node 
was described by the function of /- type [33]. Thus phenomenological approach enables to 
describe experimental nodal structure of SOP by different types of functions. It should be 
noted that making use of the reduction of basis functions of rotational group onto the group 
D 2 h, linear combinations of odd spherical functions can be related to the functions of Table 
4. Whereas, since we used general rules for construction of two-electron wavef unctions, our 
conclusions on te nodal structure of each IR don't depend on the basis functions. 

IV. CONCLUSION 

Triplet two-electron wavefunctions for the space groups with point groups Th and D 2 h are 
constructed. It is shown that total basis for the group includes each one-dimensional IR three 
times and each three- dimensional IR nine times and the basis functions of each IR vanish 
at different directions. Thus in general case there are no symmetry reasons for vanishing 
of SOP of any particular symmetry at any direction (Blount theorem [13]). Spontaneous 
magnetic fields (time-reversal symmetry violations) result in reduction of symmetry group 
of spin part to D 2 h- In this case there are two physically non-equivalent basis sets with zero 
projection of total spin on quantization axis : basis of IR Au wit four nodes in directions 
±X and ±Y, basis B lu + B 2u with two nodes in directions ±Z and basis B 3u with nodes 
in directions ±X, ±Y and ±Z. Basis A u with four nodes (or basis B 3u with six nodes [8]) 



A B (k) = A(l - kf) 



(11) 
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corresponds to A-phase of PrOs±Sb\2 and basis B\ u + £? 2lt with two nodes corresponds to the 
B— phase. Thus group theoretical account for the time-reversal symmetry violation made 
possible to connect experimentally observed nodal structure of SOP in PrOs^Sb^ with 
the index of IR of the symmetry group without additional assumptions on the shape of basis 
functions. 
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Table 1. Symmetry of two-electron wavefunctions for the direction k z (group T h ) 



Basis 


IR of D 2h 


IR of T h 


^ z 




Alg, Elg, E 2g 


$* 

z 


B\ u 


T 




B 2u 


T 




B 3u 


T 




A u 


Aiu, E lu , E 2u 



Table 2. Triplet two-electron wavefunctions for the star {k x , k y , k z } (group T h ) 
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IR of T h 


Function 


IR Of L>97, 




$* V 




1 u 


<f>* y 






& x 


J - / 2u 




^ x^ 




x u 


$' X 


It 1 „. 






Oil 


A lu 


$*z + $*y + 


A u 


Em 




A u 


Eiu 


$* z + £ 2 $*y + e& x x 


A u 



Table 3. Triplet two-electron wavefunctions for the group Th. 



IR 


Function ' 


nodes 


A»(£) 


+ ($* 9 +$* 11 -$ t 12 -$ t 10 )i/ + 




4i(y) 


_$*)^ + ($* 9 -$« 11 +$* 12 -$« 10 )i + (** 5 -$|+s«-**)z 


x,z 


4,0*) 












T-(x) 


($ 9 +^n+$i 2 +^io)y 




($' 5 +$^+$*+$ i 7 )5 






x,y 


($i-$* 11 -$' 12 +$' 10 )i 


($' 5 -$*-$*+$*)£ 


t:( Z ) 


($*-$ 2 +$|-^)i 


x,z 


($<-$' 11+ $' 12 -$^ 


($;-$*+$«-$«)y 




($J+$|-*«-$«)y 




($*+$*-$ t 3 -$*)z 


($ 9 +$* 11 -$* 12 -$ t 10 )x 



*) Function Q\ is denned by formula (7), the functions are the results of action of 
rotations of group Th in Kovalev [32] notations. 

Table 4. Triplet two-electron wavefunctions for group D 2 h- 
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IR Do h 

iil J - x Zfl 


Function 


Nodes 


ill 


VI 2 3 1 4/ 


x, y 


Bin 


($'-$*+$*-<l>i)z 




B2u 




yz 


B 3u 


($ < 1 + $*+$*+$*)5 


X,Y,Z 



